3-D General Relativistic MHD Simulations of Generating Jets

We have performed a first fully 3-D GRMHD simulation with Schwarzschild black hole with a free falling corona. The initial simulation results show that a jet is created as in previous axisymmetric simulations. However, the time to generate the jet is longer than in the 2-D simulations. We expect that due to the additional azimuthal dimension the dynamics of jet formation can be modified.Comment: 4 pages Proc. Oxford Radio Galaxy Workshop ed. R. Laing & K. Blundell (San Francisco: PASP) in press (revised

Acceleration Mechanics in Relativistic Shocks by the Weibel Instability

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks may be responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated long-term particle acceleration associated with relativistic electron-ion or electron-positron jet fronts propagating into an unmagnetized ambient electron-ion or electron-positron plasma. These simulations have been performed with a longer simulation system than our previous simulations in order to investigate the nonlinear stage of the Weibel instability and its particle acceleration mechanism. The current channels generated by the Weibel instability are surrounded by toroidal magnetic fields and radial electric fields. This radial electric field is quasi stationary and accelerates particles which are then deflected by the magnetic field.Comment: 17 pages, 5 figures, accepted for publication in ApJ, A full resolution ot the paper can be found at http://gammaray.nsstc.nasa.gov/~nishikawa/accmec.pd

Simulation of electrostatic ion instabilities in the presence of parallel currents and transverse electric fields

A spatially two-dimensional electrostatic PIC simulation code was used to study the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a field-aligned drift for L is much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field aligned drift and the E x B drift are individually subcritical. The simulation results show that the growing ion waves are associated with small vortices in the linear stage, which evolve to the nonlinear stage dominated by larger vortices with lower frequencies

Particle Acceleration, Magnetic Field Generation, and Emission in Relativistic Shocks

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating into an ambient plasma. We find small differences in the results for no ambient and modest ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. The small scale magnetic field structure generated by the Weibel instability is appropriate to the generation of ``jitter'' radiation from deflected electrons (positrons) as opposed to synchrotron radiation. The jitter radiation resulting from small scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks.Comment: 6 pages, 1 figure, revised and accepted for Advances in Space Research (35th COSPAR Scientific Assembly, Paris, 18-25 July 2004

Weibel instability and associated strong fields in a fully 3D simulation of a relativistic shock

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) excited in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a new 3-D relativistic particle-in-cell code, we have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. The simulation has been performed using a long simulation system in order to study the nonlinear stages of the Weibel instability, the particle acceleration mechanism, and the shock structure. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic (HD) like shock structure. In the leading shock, electron density increases by a factor of 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. We discuss the possible implication of our simulation results within the AGN and GRB context.Comment: 4 pages, 3 figures, ApJ Letters, in pres

From electrons to Janskys: Full stokes polarized radiative transfer in 3D relativistic particle-in-cell jet simulations

The underlying plasma composition of relativistic extragalactic jets remains largely unknown. Relativistic magnetohydrodynamic (RMHD) models are able to reproduce many of the observed macroscopic features of these outflows. The nonthermal synchrotron emission detected by very long baseline interferometric (VLBI) arrays, however, is a by-product of the kinetic-scale physics occurring within the jet, physics that is not modeled directly in most RMHD codes. This paper attempts to discern the radiative differences between distinct plasma compositions within relativistic jets using small-scale 3D relativistic particle-in-cell (PIC) simulations. We generate full Stokes imaging of two PIC jet simulations, one in which the jet is composed of an electron-proton ($e^{-}$-$p^{+}$) plasma (i.e., a normal plasma jet), and the other in which the jet is composed of an electron-positron ($e^{-}$-$e^{+}$) plasma (i.e., a pair plasma jet). We examined the differences in the morphology and intensity of the linear polarization (LP) and circular polarization (CP) emanating from these two jet simulations. We find that the fractional level of CP emanating from the $e^{-}$-$p^{+}$ plasma jet is orders of magnitude larger than the level emanating from an $e^{-}$-$e^{+}$ plasma jet of a similar speed and magnetic field strength. In addition, we find that the morphology of both the linearly and circularly polarized synchrotron emission is distinct between the two jet compositions. We also demonstrate the importance of slow-light interpolation and we highlight the effect that a finite light-crossing time has on the resultant polarization when ray-tracing through relativistic plasma.Comment: 21 pages, 13 figures; accepted for publication in A&

Particle Acceleration in Relativistic Jets due to Weibel Instability

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. While some Fermi acceleration may occur at the jet front, the majority of electron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation (Medvedev 2000) from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.Comment: ApJ, in press, Sept. 20, 2003 (figures with better resolution: http://gammaray.nsstc.nasa.gov/~nishikawa/apjweib.pdf

Particle Acceleration and Magnetic Field Generation in Electron-Positron Relativistic Shocks

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-positron jet front propagating into an ambient electron-positron plasma with and without initial magnetic fields. We find small differences in the results for no ambient and modest ambient magnetic fields. New simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. Furthermore, the non-linear fluctuation amplitudes of densities, currents, electric, and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper at the comparable simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. Additionally, we have performed simulations with different electron skin depths. We find that growth times scale inversely with the plasma frequency, and the sizes of structures created by the Weibel instability scale proportional to the electron skin depth. This is the expected result and indicates that the simulations have sufficient grid resolution. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform, small-scale magnetic fields which contribute to the electron's (positron's) transverse deflection behind the jet head.Comment: 18 pages, 8 figures, revised and accepted for ApJ, A full resolution of the paper can be found at http://gammaray.nsstc.nasa.gov/~nishikawa/apjep1.pd

Particle Acceleration and Radiation associated with Magnetic Field Generation from Relativistic Collisionless Shocks

Shock acceleration is an ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet front propagating through an ambient plasma with and without initial magnetic fields. We find only small differences in the results between no ambient and weak ambient magnetic fields. Simulations show that the Weibel instability created in the collisionless shock front accelerates particles perpendicular and parallel to the jet propagation direction. The simulation results show that this instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.Comment: 4 pages, 1 figure, submitted to Proceedings of 2003 Gamma Ray Burst Conferenc

Particle Acceleration, Magnetic Field Generation, and Associated Emission in Collisionless Relativistic Jets

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., active galactic nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs within the downstream jet. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons has different properties than synchrotron radiation which assumes a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.Comment: 4 pages, 3 figures, contributed talk at the workshop: High Energy Phenomena in Relativistic Outflows (HEPRO), Dublin, 24-28 September 2007. Fig. 3 is replaced by the correct versio